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AXEL MALQVIST AND DANIEL PETERSEIM 

Abstract. We discuss the possibility of computing eigenpairs 

of some prototypical linear second-order self-adjoint elliptic par- 
tial differential operator (or its high-resolution finite element dis- 
cretization) by numerical upscaling techniques. We compute a low- 
dimensional generalized finite element space that preserves small 
eigenvalues in a superconvergent way. The approximate eigenpairs 
are then obtained by solving the corresponding low-dimensional 
algebraic eigenvalue problem. The rigorous error bounds are based 
on two-scale decompositions of (fl) by means of a certain Clement- 
type quasi-interpolation operator. 



1. Introduction 

In this paper, we present and analyze a numerical upscaling tech- 
nique for computing eigenpairs of self- adjoint linear elliptic second or- 
der differential operators with arbitrary positive bounded coefficients. 
The precise setting of the paper is as follows. Let Q c be a 
bounded Lipschitz domain with piecewise flat boundary and let A e 
L°° {^l, Mfyj^) be a matrix- valued coefficient with uniform spectral bounds 
< a < /3 < oo, 

(1.1) a{A{x)) C [a, 13] for almost all a; e 

Consider the prototypical self-adjoint second order linear partial differ- 
ential operator A : IIo(^) — > H^^{Q) defined by 

(1.2) Av := - dw{AVv) for v E HJ(1]) 

and the associated eigenproblem: find pairs consisting of an eigenvalue 
A e R and associated non-trivial eigenfunction u E V := H^{Q) such 
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that 

(1.3) Au = Xu. 

We are mainly interested in the small eigenvalues of A or some ac- 
curate approximation Ah- The discretized operator Ah shall be related 
to a conforming finite element space of dimension Nh based on some, 
possibly very fine, mesh Th of width h. 

Popular approaches for the computation of these eigenvalues include 
Lanczos/Arnoldi-type iterations (as implemented, e.g., in |LSY98j ) or 
the QR-algorithm applied directly to the A/'/i-dimensional finite element 
matrices. If more structure can be exploited (e.g., a hierarchy of finite 
element meshes and/or spaces) some preconditioned outer iteration for 
the eigenvalue approximation may be performed and linear problems 
are solved (approximately) in every iteration step |Hac79] . |KN03bj . 
|KN03aj : see also |Ney03| and |BBS08j and references therein. 

Our aim is to avoid the application of any eigenvalue solver to the 
large-scale matrix Ah directly. Instead, inspired by |MP11] . we will 
compute some low-dimensional approximation Ah of A (resp. Ah) 
first. This preprocessing step is done by (approximately) inverting 
the operator (resp. its discretization) for special right hand sides and 
subject to certain linear constraints. Having performed Nh <^ Nh 
of these computations, the solution of a low-dimensional Nh x A'^^^ 
eigenvalue problem by standard algebraic solvers yields approximations 
of the first Nh eigenpairs. 

We emphasize that the linear problems to be solved in our method 
are completely independent of each other. They can be computed in 
parallel without any communication. Moreover, their solutions give rise 
to some low-dimensional discrete approximation space and the corre- 
sponding low-dimensional discrete operator Ah approximates the small 
(leading) eigenvalues and, hence, the macroscopic response of the op- 
erator, accurately. In this sense, our method may be understood as 
numerical upscaling or computational homogenization. 

Our method is related to some coarse finite element mesh with maxi- 
mal width H. The accuracy of the approximate eigenvalues is expressed 
in terms of H. Without any assumptions on the smoothness of eigen- 
functions, we prove that the error scales like H'^. Note that a standard 
first-order conforming finite element computation yields under full 
H^(f2) regularity, see e.g. |LarOO ]. Under such strong assumption of 
high regularity, the two-grid method of |XZ01j allows certain postpro- 
cessing (solution of linear problems on the fine scale) of the coarse finite 
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element eigenpairs to increase the accuracy to if^. However, the regu- 
larity assumption is essential and not justified on non-convex domains 
or for heterogeneous and highly variable coefficients. 

Since our estimates are both, of high order (at least H"^) and indepen- 
dent of the underlying regularity, the accuracy of our approximation 
may actually suffice to fall below the error of the fine scale discretization 
which scales algebraically in the fine mesh size h but depends on the 
regularity of the data (convexity of Q, differentiability and variability 
of A) in a crucial way. 

In cases with singular eigenfunctions (due to re-entrant corners in the 
domain or isolated jumps of the coefficient), it may be advantageous to 
use modern mesh-adaptive algorithms driven by some a posteriori error 
estimator as proposed and analyzed, e.g., in |LarOO] . |DPR03j . |CGllj . 
|GMZ09j . |GG09j . |MMllj . ^Ul2\ . We are not competing with these 
efficient algorithms. However, adaptive mesh refinement has its limi- 
tations. For instance, if the diffusion coefficient A is highly variable, 
the mesh width has to be sufficiently small to resolve the oscillations 
|PS12j . For problems in geophysics or material sciences with charac- 
teristic geometric features on microscopic length scales, this so-called 
resolution condition is often so restrictive that the initial mesh must 
be chosen very fine and further refinement exceeds computer capacity. 
Our method is especially designed for such situations which require 
coarsening rather than refinement. 

The main tools in this paper are decompositions of H^(f2) into coarse 
and fine parts which are presented in Section [2j The method for the 
approximation of eigenvalues is presented and analyzed in Section |3] 
which further discusses practical aspects and complexity issues. Nu- 
merical experiments are presented in Section |4| 

In the remaining part of this paper, we will frequently make use of 
the notation a < 6 which abbreviates a < Cb, with some multiplicative 
constant C > which only depends on the domain Q and the parameter 
p (cf. ( |2.1 ) below) that measures the quality of some underlying finite 



element mesh. We emphasize the C does not depend on the mesh size, 
the eigenvalues, and the coefficient A. Furthermore, a ^ b abbreviates 
a S.b < a. 



2. Two-scale Decompositions 

Two-scale decompositions of functions u ^ V = HQ(fi) into some 
macroscopic/coarse part Ucs G Ks plus some microscopic/fine part Wfs G 
Vis with a certain orthogonality relation are at the very heart of this 
paper. We understand macroscopic or coarse in the sense that Vcs 
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is a low-dimensional generalized finite element space based on some 
coarse finite element mesh. We understand microscopic or fine in the 
sense that Vfs contains highly oscillating functions which cannot be 
represented on the coarse mesh. The decompositions will be presented 



in Section 2.4 after some preliminaries presented in Sections 2.1 2.3 



2.1. Finite Element Mesh. We denote the underlying coarse mesh 
by Tff. The mesh is assumed to be some regular (in the sense of |Cia87] ) 
finite element mesh of fl into closed simplices with mesh-size function 
0<H eL°°{n) defined hj H\t = diam(T) =: Ht for all T e T. The 
mesh size may vary in space. However, we will not exploit the possible 
mesh adaptivity in this paper. The error bounds, typically, depend 
on the maximal mesh size | l^oo^-j^^. If no confusion seems likely, we 
will use H also to denote the maximal mesh size (instead of writing 
I l-^l '^^^ of capital letter H for the mesh-size instead 
of the standard choice h shall indicate that we consider a coarse mesh 
and that we are not interested in asymptotics as the mesh size tends 
to zero. As usual, the error analysis depends on the constant p > 
which represents the shape regularity of the finite element mesh Th', 

diamT 

(2.1) p := max px with pT := -r. ;=;- for T G Th, 

t&Th diam Bt 

where Bt denotes the largest ball contained in T. 

2.2. First-Order Conforming Finite Elements. The first-order con- 
forming finite element space corresponding to Th is given by 

(2.2) 

Vh '■= {v & V \ VT G Th, v\t is a polynomial of total degree < 1}. 

Let M denote the set of interior vertices of Th- For every vertex z G N , 
let 02 G Vff denote the corresponding nodal basis function (tent/hat 
function) determined by nodal values 

(pziz) = 1 and 02 (y) = for all y 7^ z G M. 

These nodal basis functions form a basis of Vh- The dimension of Vh 
equals the number of interior vertices, 

Nh := dimVy = \M\. 

2.3. Quasi-Interpolation. The key tool in our construction will be 
the bounded linear surjective Clement-type (quasi-)interpolation oper- 
ator Xh : V ^ Vh presented and analyzed in jCar99j . Given f G V^, 
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'^hv := X]zeA^(2^Ji"^)(-2)02 defines a (weighted) Clement interpolant 
with nodal values 

(2.3) {Xhv){z) := 

for z G A/". The nodal values are weighted averages of the function over 
nodal patches Uz := supp02. 

Recall the (local) approximation and stability properties of the in- 
terpolation operator Zh |Car99] : There exists a generic constant Cx^ 
such that for all v E V and for all T G Th it holds 

(|2^a) H^'^Wv - Xhv\\l'2(t) + W'^iv - Ihv)\\l2(t) < Cx„\\Vv\\l2{u;t), 



where ut '■= U{K E Th \ T n K (/}}. The constant Cj^ depends 
on the shape regularity parameter p of the finite element mesh Th (see 



(2.1) above) but not on Ht- Note that there exists a constant Coi > 
that only depends on p such that the number of elements covered by 
Ut is uniformly bounded (w.r.t. T) by Coi- Both constant, Cx„ and 
Col, may be hidden in the notation "<" introduced at the end of the 
Introduction on page [3j 

2.4. Decompositions. We define the finescale space 

Vfs := kernel Xh, 

which will take over the role of the microscopic/fine part in all subse- 
quent decompositions. 

Our particular choice of a quasi-interpolation operator gives rise to 
the following orthogonal decomposition. We write ||-|| := (•, •) to 
abbreviate the L^-norm on Q. 

Lemma 2.1 (L^-orthogonal two-scale decomposition). Any function 
V E V can be decomposed uniquely into the sum of f ^ := Xh\v^ {Xh{v)) E 
Vh and v^^^ := v — vh EVfs with 

(2.5) {vh,Vis)l'2{q.) = ^- 

The orthogonality implies stability in the sense of 

II ii2 II ii2 iiii2 
\\Vh\\ + t-fs = f . 



2.1 



Remark 2.1 (L^-projection onto the finite element space). Lemma 
shows that the L^-orthogonal projection Ily^ : V — )■ Vh onto the fi- 
nite element space Vh may be characterized via the modified Clement 



interpolation (2.3), 
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Proof of Lemma 2J_. It is easily verified that the restriction of Xh on 
the finite element space Vh is invertible. This yields the decomposition. 
For the proof of orthogonality, let Vh = J2zeAf'"H{z)(f)z & Vh and Vfs G 
Vfs. Since X//(ffs) = 0, we have that (02, ffs) = ZH{vis){z) = for 
all z G Af. This yields 

and shows that Vh and Vfs are orthogonal subspaces of V. □ 

as 



We may rewrite Lemma 2.1 



(2.6) V = Vh® Vs and {Vh, V^) = 0. 

The orthogonalization of this decomposition with respect to the 
scalar product a yields the definition of a modified coarse space V^g, 
that is the a-orthogonal complement of Vfs in V. Given v & V, define 
its the a-orthogonal fine scale projection operator G Vfs by 

aC^VjW) = a{y,w) for all w G Vfs. 

We define the energy norm |||-||| := ^/a{^^ (the norm induced by the 
scalar product a). 

Lemma 2.2 (a-orthogonal two-scale decomposition). Any function G 
V can be decomposed uniquely into v = Vcs + ffs, where 

Vcs := (1 - d)v G Ks := (1 - 

and 

ffs := G Vfs = kernel Xh- 
The decomposition is orthogonal 

(2.7) a(wcs,^^fs)L2(Q) = 0, 
and, hence, stable in the sense of 

(2.8) lll^cslll^ + lll^'fslll^ = lll^^lll^ • 
In other words, 

(2.9) V = Ks © V^fs and a(Vcs, V^fs) = 0. 

We aim to compute approximate eigenvectors in the space V^g. The 
orthogonalization has preserved the dimension of the coarse part of the 
decomposition, i.e., 

dim Vcs = dim Vh = Nh 
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Moreover, the images of the nodal basis functions under the a-orthogonal 
projection (1 — 5) onto Vcs yield a basis of V^s, 

(2.10) V,, = span{{l-^)(p.\ze^^}. 

In order to use this space for the approximation of eigenvalues and 
eigenvectors of (3.1 ), we need to approximate Nh solutions ipz 
Vis of 



(2.11) 



a{ipz,v) = a{(f)z,v) for all v G Vfg. 



These problems are linear and independent. The only difference with 
a standard Poisson problem is that there is some linear constraint hid- 
den in the space Vfs, that is, the quasi-interpolation of trial and test 
functions vanishes. 

The main reason why this preprocessing of the coarse space V^s can 
be useful is the following L^-quasi-orthogonality of the a-orthogonal 
decomposition (2.9). 

Theorem 2.3 (L^-quasi-orthogonality of the a-orthogonal decomposi- 
tion). The decomposition V = Vcs (BVis from Lemma 2.2 is L"^ -quasi- 
orthogonal in the sense that for all Vcs G and all G Vfs, it holds 

(2.12) {vcs,V{s)mn) < H^W'^VcsW \\Vv{s\\ < a~^H^ \\\vcs 



\V{s 



The decomposition is stable in the sense that 



(2.13) 



\H^'vis\?< 



a 



Proof. Given any t>cs G Vcs and Vis G \4s. Lemma |2.1| the Cauchy- 
Schwarz inequality, (2. 4, a), and (2.1) yield 

(^^cs, Vfs) = (Vc 



since {IhVcs,V{s 
mate shows 



,s - IhVcs, Vis - luVis) < H'^ 1 1 Vt;cs| I 1 1 VWfsl I , 
which is the quasi-orthogonality. A similar esti- 



{Hvis, Hvis) = {H{v{s - IhV{s), H{vis - InVis)) 

TeTH 
< a-' 



|2 



This, Friedrichs' inequality 



< IT diamfi llVfr 



and (2.8) readily prove the stability estimate. 



□ 
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3. Approximation of Eigenvalues and Eigenvectors 
This Section presents a new scheme for the approximation of eigen- 



values and eigenvectors of (1.3) and its rigorous error analysis. 



3.1. Variational Formulation and Galerkin Approximation. Re- 



call that the eigenpairs of (1.3) are characterized equivalently as the 



solution to the variational problem 

(3.1) a{u^'\v) := / {AW^) ■ Vv = X^'\u'^'\v) for all v E V, 



Jn 

where (■, ■) denotes the standard inner product in L^(f2). Recall that, 
for this problem, there exists a countable number of eigenvalues. More- 
over, since A is symmetric, all eigenvalues A^^-* {£ e N) are real and 
positive, 

< A(i) < A(2) < A(3) < . . . . 

Depending on the actual domain Q and the coefficient A, there may be 
multiple eigenvalues. A multiple eigenvalue is repeated several times 
according to its multiplicity in the enumeration above. Let u^^^ {£ E N) 
be normalized to one in L^(f2), i.e., {u^^\u^^^) = 1. It is well known 
that the eigenfunctions enjoy (or, in the case of multiple eigenvalues, 
may be chosen such that they fulfill) the orthogonalities 

(3.2) a(M('),M(™)) = (m(^),m(™)) = iff^m. 



The Galerkin discretization of (3.1) with respect to the generalized 
finite element space V^s reads: find A^"* E M and non-trivial uif E Ks 
such that 

(3.3) a{u^fj,v) = AgVif, ^) for all v E Ks- 

Since Vcs is a finite-dimensional subspace of V, we can order the discrete 
eigenvalues similar as the original ones 

0<A(;)<Ag)<Ag)<---<A(f-). 

Again, multiple eigenvalues are repeated according to their multiphc- 
ity. Let also ui^J (£ = 1,2,..., Nh) be normalized to one in I/^(f2), i.e., 
(u'^sjuii^) = 1. The discrete eigenfunctions satisfy (or, in the case of 
multiple eigenvalues, can be chosen such that they satisfy) the orthog- 
onalities 

(3.4) a(«^:\«ir^) = («if,«^r^) = iii^m. 
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In the subsequent paragraphs we will present error bounds for the 
approximate eigenvalues and eigenvectors based on the variational tech- 
niques from |SF73] (which are based on |BdBSW66] on their part); see 
also IBofinj . 

3.2. Two-Scale Decomposition Revisited. The eigenfunctions al- 
low a different (with respect to Section [2]) characterization of a macro- 
scopic function, that is, any function spanned by eigenfunctions related 
to the i smallest eigenvalues. Define 

(3.5) Ei:=span{u^^\...,u'-^^}. 



We will have a closer look at the quasi-orthogonality results in Lemma [272 
given some macroscopic function u G Eg. 

Corollary 3.1 (L^-quasi-orthogonality of the a-orthogonal decompo- 
sition of macroscopic functions). Let £ E N and let u = Ucs + Wfs ^ Ei 
with \\u\\ = 1, where Ucs G Ks (resp. -Ufg G Vfs) denotes the coarse 
scale part (resp. fine scale part) of u according to the a-orthogonal 



decomposition in Lemma 2.2 Then it holds 



(3-6) lll^icslll < 

(3.7) 1 1 Kill < J£(AW)3(^) , and 




(3.8) (McMfs) < Vi 

Proof. Let 6j < 1, j = 1,2, be the coefficients in the repre- 



sentation of u by eigenfunctions, that is, u = X]j=i '^j'"''"'^- Then (3.6) 
follows from the fact that (1 — S^) is a projection and the obvious bound 



\u\ 



equation (3.1) 



For the proof of (3.7), we employ some algebraic manipulations and 



(3.9) |||Mfs|||^ = a{u,uis) = ^5ja(u(^\Mfs) = '^5j\^^\u^^\uk). 

j=i i=i 



Theorem 12.31 shows that 
(3.10) {u^'\ui,)< (^) \\\u^^^\ 



The combination of ([379l)-([3lo| together with ||k(^)|||^ = A^^^ < A^') 



and 5j < 1 yields the upper bound of |||Mfs|||. 
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The combination of Theorem 2.3 and the bounds (3.6)-(3.7) readily 
yields Q. □ 

Remark 3.1. In certain cases, e.g., if Q is convex and the coefficient 
A is constant, we have that any macroscopic function u & Ee is in 
H^(i7) and ||V^m|| < A^^^^a 1 1^1 1. Such an instance of regularity gives 
rise to an additional power of if/A*-^^ in the estimates (3.7) and (3.8) 
in Corollary 3.1 This is due to the possible modification of (3.10), 

(w^^^Mfs) = (m^^^ -Inu'-^Kufs -X^Mfs) < — lllMfslll • 



a2 



3.3. Estimates for Approximate Eigenvalues. We first introduce 
the Rayleigh quotient, which is defined for non-trivial v E V hy 

Riv) := 

Recall that the ith eigenvalue may be characterized via the minmax- 
principle (which goes back to Poincare |Poi90] ) 

(3.11) AW = min max R(v), 

seSeiV)ves\{o} 

where S^iV) denotes the set of ^-dimensional subspaces of V. This 
principle applies equally well to the discrete problem (3.3), i.e., 

(3.12) X^f) = min max R(v) 

characterizes the ith discrete eigenvalue {i < Nh)- The conformity 
Vcs C V yields monotonicity 

(3.13) X^^^<X'§ for all i= 1,2,..., Nh. 

The following theorem gives an estimate in the opposite direction. 

Theorem 3.2 (Bound for the eigenvalue error). Let H he sufficiently 
small so that H < i~^^^^/^ ■ Then it holds 

aW-aw . /aw' ' 



(3.14) ^^(,) ^ I — I for all e = 1,2,..., Nh. 



a 



Proof. Recall the definition of in (3.5) and define 

:= max Km^, m^) + 2(ucs, Mfs)|, 

uGEi:(u,u)=1 

where Ucs G Ks (resp. Ufs G Vfs) denotes the coarse scale part (resp. 
fine scale part) of u G E^ according to the a-orthogonal decomposition 
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in Lemma 2.2 Theorem 2.3 and Corollary 3.1 yield 



where we have used that, under the condition H < i ^^^\/^^ ^^e 
term Kwfs, Mfs)] is of higher order when compared with |(Mcs,Mfs)|, 

|(«fs,nf.)|<£(AW)3(if/v/S)^ 

If H is chosen small enough so that cr^'' < ^, then Lemma 6.1 in |SF73] 
shows that 

A2)<(l-ag))-iAW<(l + 2a2))AW. 
Inserting our estimate for aj^'' gives the assertion. □ 

Remark 3.2 (Smallness of H). Since all our estimates are explicit with 
respect to the eigenvalues, the condition H < ^"^'^'^a/^j iii Theo- 
rem 3.2 is a qualitative condition for the smallness of H. Note that, 
typically, a condition H < 1/A'-^^ (or even more restrictive) arises 
|BBS08t ISaulOj . The relaxed condition is possible because no regu- 
larity of the eigenfunctions has been used to establish the error bound. 

Remark 3.3 (Improved eigenvalue error bound). With regard to Re- 
mark 3J^, the error bound in Theorem |3.2| may be improved provided 
that the first i eigenfunctions are regular in the sense of 1 1 V^m^-'-' | | < 
X^^^/a. The improved bound reads 

5 



(3.15) < ^m)/a ( j for alU = 1, 2, . . . , iV^. 

This improved bound might still be pessimistic in the sense that the 
error in the £th eigenvalue/vector depends on the regularity of all pre- 
vious previous eigenfunctions. The recent theory |KO06] shows that 
this is not necessarily true. Moreover, there might be smoothness also 
in the single summands of the two-scale decomposition which is not 
exploited. 

3.4. Estimates for Approximate Eigenvectors. We turn to the 
error in the approximate eigenfunctions. Again, we follow the receipt 
provided in |SF73] . 

Theorem 3.3 (Estimates of Eigenvector error). Let X^^^ be an eigen- 
value of multiplicity r, i.e. A*^^-* = ... = with corresponding 
eigenspace spanned by the orthonormal basis Le t the pair 
{A^\ tics'*} be a Galerkin approximation solving equation (3.3) such that 
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||'Wcs^||l2(q) = 1. Then there exists an eigenfunction in span{{u^^^^'>}'^^Q) , 
let us denote it u^^\ such that 

(3.16) 1 1 - I If < (AW)^£V2^ + + p)\X^^))H^^ 

Ml Ml Q,Z Q,d 

where p = maxj^{^,^+i,...,£+,_i} 

Proof. For any v G spand-u^^+^^j^'g) fulfilling ||f ||z,2(f^) = 1 it holds, 

= AW-2AW(.,«W) + a2) 

= AW(2-2(.,.W)) + a2)-AW 

(3.17) =A(^)||.-nW||i,(,) + A(J)-AW. 

It remains to bound the norm. By applying the analysis presented 
in Lemma 6.4 and Theorem 6.2 in |SF73j we have that there is a normal 
eigenfunction u^^^ E span({M'^^"'"*)}[~Q) such that, 

h^'^ -^^f < 4(1 + p)^||^fs|li.(n) < 4(i + p)^if^«-MlKs^ 

where Wfs denotes the fine scale part of u^^'^ according to the a-orthogonal 



deco mposi tion in Lemma 2.2 The theorem follows by combining equa- 
tion (3.17) with V = u^^\ Corollary 3.1, and Theorem 



3.2 



□ 



3.5. Approximation of the Ideal Coarse Space. The previous re- 
sults have theoretical impact, but the method is not yet feasible be- 
cause the definition of the coarse space Vcs involves the inversion of 
some infinite dimensional operator. That is why we consider the case 
where the original eigenvalue problem has been discretized first by some 
suitable finite element space related to some fine triangulation Th', cf. 



Section [3.5.1[ Moreover, we discuss complexity issues in Sections [3. 5. 2 
13X31 

3.5.1. Eigenvalues of High-Resolution Finite Element Discretization. 
Note that the previous estimates apply verbatim to the case, where 
V is replaced with some fine scale finite element discretization space 
Vh D Vh- One might think of Vh being the space of conforming piece- 
wise affines (or higher-order polynomials) with respect to some fine 
triangulation Th of mesh-size h < H. 

In this regard, our approach is a method for approximating the small- 
est eigenvalues of the corresponding finite element discretization by first 
solving a few large scale linear equations and then a small scale eigen- 
value problem on top of that. The accuracy of this approach is only 
linked with the coarse mesh-size H. Since our estimates are of high 
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order (at least proportional to if^) and independent of the underly- 
ing regularity, the accuracy may actually suffice to fall below the error 
of the fine scale discretization which scales algebraically in h but de- 
pends on the regularity of the data (convexity of f2, differentiability 
and variability of A) in a crucial way. 

3.5.2. Localization of Microscopic Computations and Macroscopic Com- 
pression. The construction of the coarse space V^s is based on fine scale 
equations formulated on the whole domain VL which makes them ex- 
pensive to compute. However, in |MP11] it was shown that "^ipz decays 
exponentially fast outside of the support of the coarse basis function 
02- We specify this feature as follows. Let /c G N. We define nodal 
patches u}z,k of k coarse grid layers centered around the node z G A/" by 

^z,i ■= supp (pz = U {T E Th \ z E T} , 

^^■^^^ C0z,k ■.= U{T eThITH u,,k-i ^ 0} for k>2. 

The result in the decay of ^(pz in |MP11] can be expressed as follows. 
For all vertices z E M and for all /c G N, it holds 

(3.19) \\A'''^U^l.,^,^^^^,<e-^-'^^^^^^ 

For moderate contrast this motivates the truncation of the com- 
putations of the basis functions to local patches ujz k- We approximate 



-ipz = d4>z e Vk from ( |2.11[ ) with ipz,k e Vis{uJz,k) := {v E Vfs \ vjov, 
0} such that 

(3.20) a{i)z,v) = a{(f)z,v) for all G Vfs(a;;,,fc). 
This yields a modified coarse space V^g with a local basis 

(3.21) = span{0, - tP^^k \zeM}. 

The number of non-zero entries of the corresponding stiffness matrix is 
proportional to k'^Nfj (note that we expect Nfj non-zero entries without 
the truncation). Due to the exponential decay, the very weak condition 
k ^ log H implies that the perturbation of the ideal method due to this 



truncation is of higher order and the estimates in Theorem 3^ and 3^ 
remain valid. We refer to |MP11] for details and proofs. 

3.5.3. Complexity. Finally, we shall comment on the overall complex- 
ity of our approach. Consider quasi-uniform meshes of size H resp. h 
and corresponding conforming first-order finite element space Vh and 
Vh- We want to approximate the eigenvalues related to Vh. Therefore, 
we need to solve A'^^^ linear problems with approximately k'^Nh/Nn 
degrees of freedom each; the parameter k being the truncation pa- 
rameter as above. Since almost linear complexity is possible (using. 
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e.g., multilevel preconditioning techniques), the cost for solving one of 
these problems up to a given accuracy is proportional to the number 
of degrees of freedom NiJNh up to possible logarithmic factors . This 
yields an overall complexity of k'^Nh\og{Nh) (resp. NnNhloglNh) if 
k'^ > Nh) for setting up the coarse problem. Note that this effort can 
reduced drastically either by considering the independence of the linear 
problems in terms of parallelism or by exploiting a possible periodicity 
in the problem and the mesh. In the latter case, only very few of the 
problems have to be computed because all the other ones are equivalent 
up to change of coordinates. 

On top of the assembling, an A'^f/- dimensional eigenvalue problem is 
to be solved. The complexity of this depends only on Nh, the number 
of eigenvalues of interest, and the truncation parameter k but not on 
the critically large parameter N^. 



4. Numerical Experiments 

Two numerical experiments shall illustrate our results. We focus on 
the case without localization. The localization (as discussed in Sec- 
tion 3.5.2 has been studied extensively for the linear problem |MPllj . 
In the present context of eigenvalue approximation, we are interested 
to observe the enormous convergence rate which is 4 or higher for the 
eigenvalues. In order to achieve this rate also with truncation, patches 
have to be large (at least 4 layers of elements) which pays off only 
asymptotically when H is small enough. 

4.1. Constant coefficient on L-shaped domain. Let n := (—1, 1)^\ 
[0, 1]^ be the L-shaped domain. Consider the constant scalar coefficient 
Ai = 1. Consider uniform coarse meshes with maximal mesh widths 
y/2H = 2-1, . . . , 2-^ of as depicted in Figure |ll 

The reference mesh Th has maximal mesh width h = 2~'^/\/2. We 
consider some PI conforming finite element approximation of the eigen- 
values on the reference mesh Th and compare these discrete eigenvalues 
Xf^ with coarse scale approximations depending on the coarse mesh size 
H. 

Table [T] shows results for the case without truncation, i.e., all linear 
problems have been solved on the whole of fl. For fixed i, the rate of 
convergence of the eigenvalue error A^^ — X^/^^ in terms of H observed 
in Table [T] is between 6 and 7 which is even better than predicted in 
Theorem 13.21 and in Remark 13.11 



4.2. Rough coefficient with multiscale features. Let := (0, 1)^ 
be the unit square. The scalar coefficient A2 (see Figure |2]) is piecewise 
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Figure 1. Initial uniform triangulation of the L-shape 
domain (5 degrees of freedom). 

constant with respect to the same uniform Cartesian grid of width 
2~^. Its values are taken from the data of the SPEIO benchmark, see 
http : //www . spe . org/web/ csp/. The coefficient is highly varying and 
strongly heterogeneous. The contrast for A2 is large, f3{A2)/a{Ai) ^ 
4-10^. Consider uniform coarse meshes of size \/2H = 2^^, 2^^, . . . , 
of Q as depicted in Figure [2] Note that none of these meshes resolves 
the rough coefficient A2 appropriately. Hence, (local) regularity cannot 
be exploited on coarse meshes. 

Again, the reference mesh Th has width h = 2^ \/2 and we compare 
the discrete eigenvalues A^^^^ (with respect to some PI conforming finite 
element approximation of the eigenvalues on the reference mesh Th) 
with coarse scale approximations depending on the coarse mesh size 
H. Table [2] shows the errors and allows us to estimate the average 
rate around 4 which matches our expectation from the theory. We 
emphasize that the large contrast does not affect the accuracy of our 
method in approximating the eigenvalues X^i^\ However, the accuracy 

of A)j is affected by the high contrast and the lack of regularity caused 
by the coefficient. 
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merical experiment and initial uniform triangulation of 
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